Number— phase uncertainty relations: verification by homodyning 



T. Opatrny 1 * 2 , M. Dakna 1 , and D.-G. Welsch 1 
1 Friedrich-Schiller-Universitat Jena, Theoretisch-Physikalisches Institut 
Max-Wien-Platz 1, 07743 Jena, Germany 
2 Department of Theoretical Physics, Palacky University 
Svobody 26, 77146 Olomouc, Czech Republic 
(April 18, 1997) 



: 

on . 
on ■ 

5h ! 
Oh- 

<: 
oo : 



> : 
. 

m ■ 
o ■ 
^t- : 
o . 
■ 

ON 1 

^: 

9 L,: 
^— > ■ 

G ■ 

3 : 
cr. 



It is shown that fundamental uncertainty relations between 
photon number and canonical phase of a single-mode optical 
field can be verified by means of balanced homodyne mea- 
surement. All the relevant quantities can be sampled directly 
from the measured phase-dependent quadrature distribution. 
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Although the problem of number-phase uncertainty 
has widely been studied, there has been no direct ex- 
perimental verification of fundamental uncertainty rela- 
tions (URs). What can be the best way of doing that? A 
powerful and perhaps ultimate method for measuring the 
quantum statistics of traveling optical fields has been bal- 
anced homodyne detection. The quantity that is directly 
measured is the probability distribution p(x, <&) of the 
phase-dependent quadrature x(d) = 2~ 1 / 2 (ae~ t ' a +0) e 1 ^), 
where a (&') is the bosonic annihilation (creation) opera- 
tor of the (single-mode) signal field and # corresponds to 
the local-oscillator phase. It has been shown that p(x, $) 
for all phases •& in a n interval contains all knowable in- 
formation about the quantum state of the signal field and 
can be used to reconstruct the Wigner function applying 
inverse Radon transformation techniques Q . 

Since the Wigner function is a full description of the 
quantum state, it can be used to calculate other impor- 
tant features of the field, such as the photon-number and 
phase statistics and their associated URs ||. The in- 
verse Radon transform requires a three-fold integration 
of the measured data, and the calculation of the den- 
sity matrix in the Fock basis can then be accomplished 
with two integrals. One summation eventually yields the 
photon-number moments, and one sum and one integral 
must be performed to obtain the (canonical) phase mo- 
ments. Hence, six- and seven-fold transformations of the 
recorded data are required for UR verification at least. Of 
course, large amount of data manipulation accumulates 
various errors and the physical nature of the uncertainties 
becomes less transparent. 

Recent progress has offered new possibilities of deter- 
mining the photon-number statistics in a more direct way 
avoiding the detour via the Wigner function. It has been 
shown that both the density-matrix elements g nn i in the 
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Fock basis || and the moments and correlations {a lk a k ) 
j| can directly be sampled from the recorded data ac- 
cording to a two-fold integral transform 



Here, A is the quantity that is to be determined and 
Ka{x,$) is the corresponding integral kernel (sampling 
function). In particular, A can be identified with p n = 
Qnn or (il k ) = ((a^a) k ) in order to determine the photon- 
number distribution or the photon-number moments. 

In contrast to the photon number, the phase has re- 
mained a troublesome variable. Recently it has been 
suggested to measure the canonical phase distribution 
asymptotically by replacing the local-oscillator with a 
reference mode prepared in so-called reciprocal binomial 
states - a method that is state dependent and hardly re- 
alizable at present g. It has also been suggested to use 
balanced homodyning for asymptotically sampling phase 
distributions (6). Again, the method is state dependent 
and the exact phase distributions can be obtained only 
in the sense of a limiting process. Moreover, with in- 
creasing (photon-number) excitation of the state under 
consideration an increasing number of highly oscillating 
functions must be summed up to obtain the correspond- 
ing integral kernels and high-accuracy sampling of the 
phase distributions from the recorded data becomes very 
difficult. The problem is similar to that of a calculation 
of phase quantities from the sampled density-matrix ele- 
ments in the Fock basis. With increasing excitation of the 
state a larger number of density-matrix elements must be 
sampled, so that the accumulated error may eventually 
dominate the result. 

In what follows we show that fundamental number- 
phase URs that are based on exponential phase mea- 
sures can be verified directly from the homodyne data 
according to Eq. ([!]), without making a detour via other 
quantities and without making any assumptions and ap- 
proximations with regard to the state. Let us consider 
a phase distribution p(ip) and its exponential phase mo- 
ments defined by 
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In particular for the canonical phase the relation = 
(E k ) is valid, where the exponential phase operator E is 
given by E — (n + l) _1 / 2 d (for the canonical phase and 
the exponential phase operator, see, e.g., Qj). In order 
to define a mean phase (p that is independent of the cho- 
sen phase window, the first-order exponential phase mo- 
ment has been introduced in the definition, (p = arg^i = 



1 



arg(-E'), and it has been used to define phase uncertainty 
measures, such as 



Aip = arccos |^i | = arccos | (E) | 



(3) 



[|| , the Bandilla-Paul dispersion (jg P |)| , and the Holevo 
dispersion a\ 

ctbp = sinA^j, (7h = tan Aip. (4) 

It can then be proven that the UR 

An tan Aip > | (5) 

is valid || , which is equivalent to the Holevo UR jl0| 

{Anfol > i (6) 

[(An) 2 = (n 2 ) — (n) 2 ]. Although these URs are exact, they 
are, in a sense, weak. This means that they also allow 
for such values of An and Aip for which no state exists 
(for more specific URs and the corresponding minimizing 
states, see |y],[l2]] ) . Further, URs that are based on the 
Susskind-Glogower trigonometric operators C = ^(E + 
£t) and S= ^-(E- have also been studied @,||||, 

AnAC > h\(S)\, AnAS > h\{C}\, 



AS AC > iftjo 



(7) 
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[{AC) 2 = (C 2 ) - (C) 2 , and (AS*) 2 accordingly]. Note 
that the squares of the trigonometric operators can be 
written as 



c 2 = hi 
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# 2 )-±|0)<0|, 
B+ a )-||0)(Q|. 
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From the definitions of Aip, AC and AS we see that for 
measuring them it is sufficient to measure the exponential 
phase moments "Ji = (E) and \&2 = (E 2 ) and the density- 
matrix element goo — (|0)(0|). The determination of the 
photon-number uncertainty requires measurement of (n) 
and (n 2 ). 

It is well known that goo, (n), and (n 2 ) can directly be 
sampled from the data recorded in balanced homodyn- 
ing applying Eq. (|J). The kernel Kooix,^) for goo can 
be taken from Koo(x,$) = tt -1 \, — x 2 ), where 
<i>(a, b, z) is the confluent hypergeometric function. The 
kernels K n (x, , d) and K n 2(x,$) for (n) and (n 2 ), respec- 
tively, can simply be obtained from the sampling formula 
for the normally ordered moments and correlations of 
bosonic operators 0, 
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from which we find that 
K n 2 (x, 1}) - , 
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Note that for large |a;| the leading terms in these expres- 
sions determine the kernels for the energy moments in 
the classical limit. 

Let us now turn to the problem of direct sampling of 
the exponential phase moments of the canonical phase. 
Provided that the corresponding integral kernels exist, 
their asymptotic behavior for large \x\ can be obtained 
from considering the classical limit. Since in classical 
physics the phase probability distribution p(ip) can be 
obtained from the phase-space probability distribution 
W(r, ip) according to p(ip) — J rdr W(r, ip), the expo- 
nential phase moments Vffc can be given by 



* fc = / / rdrW(r,tp)e tkv . 



(14) 



Further, the quadrature probability p(x, ■&) is given by 
the Radon transform 



p{x,-d) = I dip I rdr W(r, ip) S[x — r cos($ — ip)] . (15) 

J2-n Jo 

Let us now assume that $ can be related to p(x, $) ac- 
cording to Eq. (0) [A-r^k and K A {x,$) -> K k (x,&)]. 
Substituting in this equation for p(x, $) the result of 
Eq. Jl5"| ) and comparing with Eq. (|l4|), we observe that 
Kk{x,d) can be written as Kk(x,i)) — e lki} Kk(x), where 
Kk(x) must satisfy the integral equation 



dtp e ikv K k (r cos tp) = 1 



(16) 
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for any r > 0. From Eq. ( [lq ) we can see that Kk(x) is not 
uniquely defined. First, any function of parity (— 1) +1 
can be added to Kk{x) without changing the integral. 
Second, any polynomial of a degree less than k can also 
be added to Kk(x). As can be verified, Eq. ( fl6| ) is solved 
using the functions 



n+1 (x) = ±(-l) m (2m + 1) sign(x) (17) 



and 



K 2m (x) =7r- i (-ir+ i mlog|x| 



(18) 



(to = 0, 1, 2, . . .). It is worth noting that since Eq. jl5| ) 
is also valid when W(r, ip) is the quantum-mechanical 
Wigner function, using in homodyne detection the ker- 
nels (|l7|) and ( p8[ ) over the whole x axis would yield 
the exponential phase moments of the phase quasi- 
probability distribution defined by the radially integrated 
Wigner function. 

With regard to the canonical phase, the kernels (|l7j ) 
and (Jlq) are of course valid only for large |sc|. To obtain 



2 



them for arbitrary |zc|, we recall that in quantum physics 
\ffc can be given by 



*fe = (E k ) = ^ Qn+kn 
n=0 



(19) 



in place of Eq. ([14]). Expressing p(x,i9) in terms of the 
density-matrix elements g nn i as 



(20) 



and assuming that Eq. (|l|) applies to ^ k , we again find 
that K k (x, i9) = e M K k (x), but in place of Eq. @ iC fc (x) 
must now satisfy the integral equation 



2tt 



dx K k (x)ip n (x)ip n+k (x) = 1 



(21) 



for every n = 0,1,2,... . Here, ip n (x) = (2"n!V^)~ 1/2 
cxp(— x 2 /2)H n (x) are the energy eigenfunctions of a har- 
monic oscillator, with H n (x) being the Hcrmite polyno- 
mials. From Eq. (|21| ) and the properties of the Hermite 
polynomials |l4 1 the same arbitrariness in the determina- 
tion of K k (x) as in the classical limit is found. To derive 
an explicit expression, we use the expansion 

" ^ ■ y/n\(n + k)\ 



oo oo 

= EE 



(-1)' 



^ ^ x /n!(n + fc)! m! 



fn+m ~ n+m+fc 



(22) 



(: : introduces normal order) and apply Eq. (|ll| ) in order 
to represent in the form of Eq. ([j]). Provided that 
all the correlations (a {n a n+k ) exist, we derive 



K k (x) = (2^)- 1 ^C ; (fe) H 2;+fc (x), 



1=0 



where 



(fc) 



(Z + jfe)! 



E 



2 /+( fe /2) (2 ; + fc)! ^ Wv/( n +l)...( n +fc) 



(-1) 



(23) 
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In particular, the kernels for sampling the first two mo- 
ments can be rewritten as, on using standard relations 



KAx) = ir- 3/2 x / — o-$(2, |,-x 2 tanht) (25) 

7o Vtcosh 2 i V 2 ; V ' 

(also see jl5]]), and (apart from an irrelevant constant) 

K 2 (x) = — [°°dt ~ Io(t) — $(2, i -x 2 tanhi) (26) 
2tt 7 cosh 2 tsinht 12 J J 



(Io(t) is the modified Bessel function). Finally, it can 
be proven that the kernels exist and satisfy the condi- 
tion (pl|), i.e., we have found solutions even when the 
assumption of finite moments fails and Eq. ( |ll| ) cannot 
be used. 

The kernels K±(x) and K 2 (x) can be evaluated nu- 
merically using standard routines. They are plotted in 
Fig. |l|(a) and (b). As expected, they rapidly approach 
the classical limits given in Eqs. (0) and @ as \x\ in- 
creases. Since they differ from the classical limits only 
in a small region (of a few "vacuum-fluctuation widths" ) 
around zero, in practice their evaluation needs applying 
Eqs. fl25| ) and (^6|) [or Eq. (p4j)] only for small values 
of \x\, whereas for greater values the expressions given 
m Eqs. @ and (U) 

can be used. Note that for small 
values of \x\ power series expansion of the confluent hy- 
pergeometric function in Eqs. (|25| ) and ( ^6|) can be used. 

To verify URs connected with Aip [e.g., the relations 
(||) and (^)], the first moment of E must be measured, 
which can be accomplished with the kernel Ki(x,d). 
With regard to URs of the type given in Eqs. (0) and 
(||) , one also needs the second moment of E and the vac- 
uum probability goo- From the above, the kernels for 
sampling (C 2 ) and (S 2 ) read as 

K±{x, 0) = i [1 - $(1, i, -x 2 )] ± i cos(2d)K 2 (x) (27) 

[see Fig. 0(c)] , where K+ (x, 0) = K C 2 (x, ■&) and K_ (x, ■&) 
= K s *(x,#). 

To conclude, we have presented a method for verifica- 
tion of number-phase URs. It is based on the possibility 
of direct sampling of exponential phase moments of the 
canonical phase of a single-mode quantum state, ^ k = 
(E ), from the data recorded in balanced homodyning. 
We have shown that the corresponding kernels K k (x,i)) 
are state independent and well behaved. With increasing 
|x| they rapidly approach the classical limits. Since the 
method does not only apply to the determination of low- 
order moments, it may also be used for reconstructing the 
canonical phase distribution p(ip) as a whole. However, 
a direct (state independent) sampling of p(ip) according 
to Eq. (Q) seems to be impossible. If there would exist 
a corresponding kernel, its Fourier components with re- 
spect to 1? would be equal to the kernels for determining 
^fc- Since for chosen x the absolute values of these ker- 
nels increase with k [cf. Eqs. (Jl?]) and (|18|)1, they cannot 
be treated as Fourier coefficients of a well-behaved func- 
tion of $. Of course, this does not exclude an indirect 
reconstruction of p{ip). Measuring a limited number of 
^fc, one can use, e.g., the maximum entropy principle 
p6[ to obtain a p{ip) which best fits the measured values 
without introducing any arbitrary bias. This also offers 
the possibility of verification of URs that are not based 
on exponential phase moments |T7j , |T8| . 

This work was supported by the Deutsche Forschungs- 
gemeinschaft. 
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FIG. 1. The kernel for sampling the first (a) and second 
(b) exponential-phase moment is shown (full line) and com- 
pared with the classical result (dashed line). The kernel for 
sampling (C 2 ) (c) is shown for i9 = (1), $ = 7r/4 (2), and # 
= tt/2 (3). 
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